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The so-called ‘raspberry’ model refers to the hybrid lattice-Boltzmann (LB) and Langevin molec¬ 
ular dynamics scheme for simulating the dynamics of suspensions of colloidal particles, originally 
developed by [V. Lobaskin and B. Diinweg, New J. Phys. 6, 54 (2004)], wherein discrete surface 
points are used to achieve fluid-particle coupling. In this paper, we present a follow up to our study 
of the effectiveness of the raspberry model in reproducing hydrodynamic interactions in the Stokes 
regime for spheres arranged in a simple-cubic crystal [L. Fischer, et al, J. Chem. Phys. 143, 084107 
(2015)]. Here, we consider the accuracy with which the raspberry model is able to reproduce such 
interactions for particles confined between two parallel plates. To this end, we compare our LB sim¬ 
ulation results to established theoretical expressions and finite-element calculations. We show that 
there is a discrepancy between the translational and rotational mobility when only surface coupling 
points are used, as also found in Part I of our joint publication. We demonstrate that adding internal 
coupling points to the raspberry, can be used to correct said discrepancy in confining geometries as 
well. Finally, we show that the raspberry model accurately reproduces hydrodynamic interactions 
between a spherical colloid and planar walls up to roughly one LB lattice spacing. 


I. INTRODUCTION 

Here, we continue our investigation into the so- 
called ‘raspberry’ model, which was first introduced by 
Lobaskin and Diinweg. [T| The raspberry model allows for 
the incorporation of a colloidal particle into the lattice- 
Boltzmann (LB) algorithm by discretizing the surface 
into point particles and utilizing the coupling scheme of 
Ahlriclis and Diinweg [2] for these points. The model de¬ 
rives its name from the discretized nature of the colloid’s 
surface, which resembles a raspberry, when represented 
by molecular-dynamics (MD) beads, see Fig. [l] (right) for 
the original ‘hollow’ variant [1] and Fig. |T] (left) for the 
improved ‘filled’ raspberry introduced in Ref. |3J. 

In Ref. [S] we demonstrated the raspberry’s ability 
to accurately capture hydrodynamic interactions in a 
simple-cubic arrangement. We discussed how the point 
coupling leads to an effective hydrodynamic radius that 
can be fitted for. As shown by Ollila et al. 010 there 
is a mismatch between the effective hydrodynamic radii 
when the coupling points are placed on a spherical shell. 
Ollila et al. contend that this is related to the porosity 
of the raspberry particle, also see Part I [3], that is in¬ 
herent to this coupling method. We showed in Ref. [3] 
that one can effectively match the raspberry to results for 
solid particles within a reasonable parameter regime, via 
a procedure which we dubbed ‘filling + fitting.’ That is, 
a raspberry particle with internal coupling points (filled), 
for which effective hydrodynamic radii are established via 
hydrodynamic experiments (fitting), shows excellent nu¬ 
merical correspondence for the two fitted hydrodynamic 
radii. Moreover, the results agree well with theoretical 
predictions for a solid particle in the Stokes’ limit, in the 
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FIG. 1. (color online) Representation of the structure of the 
raspberry models used in our simulations, filled (left) and 
hollow (right), respectively. The central bead to which all 
other beads are connected via rigid bonds is shown using a 
green sphere. The blue spheres represent the beads that form 
the filled raspberry and the red ones give the surface beads 
used for the hollow variant. The radius of the beads is chosen 
to be smaller than the typical effective hydrodynamic radius 
to help visualize the internal structure. 


intermediate and long-time regime. 

In this paper, we study the confining geometry of two 
parallel plates, to further demonstrate the quality of the 
‘filling + fitting’ procedure. This geometry is of particu¬ 
lar interest, since to the best of our knowledge, no system¬ 
atic study of the quality of the particle-coupling method 
of Ref. [([] has been attempted in which the level of con¬ 
finement was the crucial parameter. A two plate geome¬ 
try has been applied in various works, see, for example, 
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Refs. [51 jS], but always with the intention to extrapo¬ 
late the results to infinite plate separations or extract 
bulk properties. For our study of small plate separations 
(strong confinement) in this manuscript, we use a variety 
of literature results to compare our results against. 

The first analytic expressions for the mobility of the 
sphere-plate system were derived by Faxen [ziiE], who 
restricted his study to spheres constrained to the mid¬ 
plane. Brenner derived an expression for mobility of a 
sphere close to a single wall [9] (in the direction per¬ 
pendicular to the wall), which was later extended to a 
parallel-plate geometry in Refs. maun. These publica¬ 
tions (EHII]) are but the tip of the iceberg when it comes 
to the semi-analytic and numerical work that has been 
carried out on these and comparable geometries. Other 
works include: semi-analytic calculations, HM5] numer¬ 
ical simulations that exploit analytic approximations to 
capture hydrodynamic behavior, [I6U2T] and simulations 
that employ coarse-grained hydrodynamic solvers to fully 
resolve hydrodynamic interactions. j22W25j In addition, 
numerous experimental studies into the diffusive behav¬ 
ior of spheres under confinement, P5H50] as well as other 
shapes, [3TH55] have been undertaken. However, it goes 
beyond the confines of this introduction to fully list and 
do justice to all these investigations. 

In the following, we use some of the above literature re¬ 
sults to prove that the hollow and filled raspberry repro¬ 
duce Stokesian hydrodynamic interactions for a sphere 
confined between two parallel plates. As in Part I, the 
‘filling + fitting’ procedure is shown to resolve the is¬ 
sue of the mismatch between the effective hydrodynamic 
radii that the hollow variant of the raspberry suffers from, 
to within acceptable numerical tolerance. We also show 
that hydrodynamic interactions are properly captured by 
the raspberry-LB method up to one lattice spacing away 
from either wall. This is a reasonable result, since lubri¬ 
cation corrections for this point-coupling method have 
not yet been established. In the right limits, we were 
able to demonstrate that the mean-square displacement 
also yields the proper position dependence of the trans¬ 
lational mobility in a thermalized LB fluid. Finally, we 
used our high-quality numerical data to establish empiri¬ 
cal relations for this position dependence for cases where 
literature values were not available or did not extend over 
the full range of plate separations. 


The remainder of this manuscript is structured as fol¬ 
lows. In Section [TT] we describe the aspects of our sim¬ 
ulation methods that were not covered in Part I [3] in 
detail. Section [II A| describes the various hydrodynamic 
experiments that we performed to determine the proper¬ 
ties of the raspberry model. For certain experiments we 
used finite-element calculations as the basis for compar¬ 
ison, the specifics of these are provided in Section |IIB| 
We provide a summary of the notations used through¬ 
out the text in Section IIC to aid the reader when going 
through the manuscript. In Section |III| we list our main 
results. These are discussed and related to previous stud¬ 
ies and our work in Ref. [3] in Section IV Finally, we give 


a summary, conclusions, and an outlook in Section |V| 


II. METHODS 

In this section, we outline the modeling approaches 
used to determine the hydrodynamic properties of a col¬ 
loid. We have split this into subsections detailing aspects 
of the hydrodynamic experiments performed to extract 
the mobility of the raspberry, the finite-element calcula¬ 
tions, and a reference list of the input parameters and 
measured quantities. For details on the construction of 
the raspberry model and the Molecular Dynamics (MD) 
and lattice-Boltzmann (LB) parameters employed in our 
investigation, we refer the reader to Ref. [5]. 


A. Hydrodynamic Experiments 

To assess the quality of the raspberry approximation 
in modeling the hydrodynamic properties of a colloid we 
performed several experiments. In the following, we use 
the term ‘quiescent’ to describe an un-thermalized (non¬ 
fluctuating, deterministic) LB fluid. We set up the sys¬ 
tem of a particle confined between two parallel plates as 
follows. A channel of height H was simulated using H- 1-2 
lattice points along the 2 -axis (one layer of boundary cells 
on either side). In the a;y-plane the box had an extent L 
in both directions and periodic boundary conditions were 
applied in these directions. We used finite-size scaling to 
remove any dependence on L, i.e., we varied the box size 
and fitted to the result for L f oo; typically the results 
were sufficiently converged for L > 5 H. The bounce-back 
boundary condition |2] was applied to the first and fi¬ 
nal layer of points in 2 -direction, leading to an effective 
hydrodynamic channel height of H with a no-slip sur¬ 
face. The height was confirmed by Hagen-Poiseuille flow 
measurements, by comparing to the theoretical fluid flow 
profile and fitting for the channel height. For all exper¬ 
iments the raspberry was initialized at a height 2 , mea¬ 
sured with respect to the center of the channel (2 = 0), 
see Fig. [2] 

The following experiments for a spherical particle con¬ 
fined between two parallel plates were performed, as de¬ 
tailed below. In Ref. [3] we discuss in detail the Reynolds 
number (Re), Peclet number (Pe), Schmidt number (Sc), 
as well as other dimensionfree constants at which these 
experiments were carried out. For the purposes of this 
manuscript, it suffices to know that the typical values of 
these numbers are Re < 0.1, Pe > 100, and Sc > 1000. 

• A velocity experiment in a quiescent fluid, see 
Fig- da). The particle was kept fixed at the initial 
position. Its surface velocity was held constant, by 
resetting the velocity in every time step. Thus an 
effective velocity boundary condition is obtained, 
without the particle changing position. The mo¬ 
bility was subsequently determined by measuring 
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FIG. 2. (color online) Visualization of the various hydrody¬ 
namic experiments carried out between two parallel plates, 
separated by a distance H. Here, a two-dimensional (2D) 
cross-section of the setups is shown. The blue arrows and 
symbols denote quantities applied to the fluid and raspberry, 
the red arrows and symbols indicate measured quantities, the 
green arrows indicate a spatially fixed raspberry. The black 
arrows signify a thermalized fluid. We refer to the text for 
a description of the experiments and applied and measured 
quantities. 


the force on the particle in the stationary state. 
This approach is similar to the one employed for the 
finite-element calcultions in Section III Bl Since the 
particle does not move, the position dependence of 
the translational mobility parallel and perpendicu¬ 
lar to the wall could be accurately assessed. 


• A force experiment in a quiescent fluid, see 
Fig. [2])b) . The particle was allowed to move freely 
and a constant force was applied to its center of 
mass. The position dependent mobility for move¬ 
ment parallel to the walls could be determined by 
extracting the terminal velocity. This terminal ve¬ 
locity was averaged over at least 10 periods of the 
(slight) oscillation in this parameter due to inter¬ 
polation artifacts. It proved unnecessary to fix the 
particle in the 0 direction, since there was no appre¬ 
ciable vertical drift. The presence of the boundary 
breaks the symmetry of the system, when the par¬ 
ticle is not in the center of the channel. Therefore, 
there can be rotational and translational cross- 
coupling terms in the hydrodynamic mobility ten¬ 
sor (HMT). [35) To examine this we also measured 
the angular velocity of the particle. 

• A torque experiment in a quiescent fluid, see 
Fig-! c). A constant torque was applied to the 
particle about an axis parallel (Fig. [2])c), left) or 
perpendicular (Fig. |2])c) , right) to the walls and the 
angular velocity was measured. Again it proved un¬ 
necessary to fix the particle in the z direction. For 
rotation about an axis parallel to the wall, we also 
fixed the particle’s position, to examine the effect 
of cross coupling. 

• A mean-square displacement (MSD) experiment in 
a thermalized fluid, see Fig. [2])d). The system was 
equilibrated until the particle fluctuated with the 
proper imposed fluid temperature. Then over many 
integration cycles (typically 10 8 per run) and for 
several runs, we kept track of the position of the 
particle. The positions were binned according to 
the following procedure. Each bin was made half 
a lattice spacing (0.5er) in width for the direction 
of interest. All sub-trajectories of our runs that 
originated in a specific bin were assigned to that 
bin and followed until they exited the bin. To en¬ 
sure a minimum smearing out of the position de¬ 
pendence, we selected only those parts of the tra¬ 
jectories that went through a slice O.ler in diameter 
around the center of the bin. This eliminated the 
trajectories that moved in and out of the bin but 
never crossed the center. For each selected trajec¬ 
tory, we determined the MSD and averaged this 
over all the relevant trajectories. Each individual 
MSD was obtained by examining such a trajec¬ 
tory over 50,000 integration steps from its ‘start¬ 
ing’ point. Because we used a low temperature 
(fcs T = O.Ole), we were able to access the long¬ 
time diffusion regime without the particle travel¬ 
ing beyond the confines of the relevant bin, despite 
the long MSD sampling length, see Part I. [3] This 
is necessary, as the Aldrichs and Diinweg coupling 
method (5: used to thermalize the raspberry parti¬ 
cles, does not accurately reproduce the short-time 
diffusion that is present for a system that satisfies 
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Stokes’ equations for all times, see the discussion 
in Ref. [3]. 


B. Finite-Elements Calculations 

In order to compare the quality of our results for the 
raspberry particles to reference curves, we utilized liter¬ 
ature values from analytic and numerical studies when¬ 
ever possible. However, for the sphere confined between 
two plates, literature results were often not available, 
to the best of our knowledge. Therefore, we performed 
finite-element calculations using the COMSOL 4.4 Multi- 
Physics Modeling Software to establish reference data. 

We set up a system with a spherical colloid of radius 
R = 0.5 pm. The sphere was located on the symmetry 
axis of a cylindrical fluid domain of height H , bounded on 
top and bottom by two no-slip plates (the bases of the 
cylinder). The cylindrical domain was bounded on its 
side by open boundary conditions (zero normal stress) 
for the fluid. We solved Stokes’ equations for an incom¬ 
pressible fluid with density p = 10 3 kgnU 3 and dynamic 
viscosity r) = 10 -3 kg nr -1 s _1 on this domain. 

To determine the various mobilities we imposed the 
following boundary conditions on the sphere. We solved 
for the force and torque exerted by the fluid and plates 
on the sphere, to which we applied a constant (angular) 
velocity boundary condition, to obtain the translational 
and rotational mobility, respectively. To ensure that the 
laminar flow condition is satisfied, we used low values of 
the (angular) velocity of the sphere: v = 10~ 6 ms -1 and 
u! = 10 -6 s . Typically, we found Reynolds numbers 
Re T and Re R < 10“ 5 . 

We varied the position of the sphere with respect to 
the plates along the cylindrical domain’s axis keeping 
the sphere’s position fixed in the other directions to de¬ 
termine the z dependence of the mobility. Similarly, we 
examined the plate-separation dependence H of the mo¬ 
bility by keeping the sphere fixed in the center of the 
cylinder and moving the plates outward. The diameter 
of the cylinder was varied to remove finite-size effects, 
but we found that a diameter of 10 H was sufficient to 
guarantee convergence. 

Whenever possible, we exploited the rotational sym¬ 
metry of the system to minimize the number of elements 
required, i.e., for perpendicular movement and rotation 
about an axis perpendicular to the wall. For the full 
three-dimensional (3D) geometries, constructing a suffi¬ 
ciently refined mesh that contained a manageable num¬ 
ber of elements proved difficult. We used local refinement 
with elements 1 nm in size over the surface of the particle 
and in disks with radius R on the bases of the cylinder 
around the cylinder’s axis. The size of the elements was 
allowed to expand radially outward. We verified that 
integration over the surface yielded the surface area to 
within a fractional deviation of 10 . We therefore ex¬ 

pect that similar integration of the force and torque over 
the surface would show limited deviation between the 


true solution of the boundary value problem and our nu¬ 
merical result. To further demonstrate the quality of our 
approximate solution, we performed mesh refinement for 
a selected number of configurations and found no signif¬ 
icant change. 

It should be noted that there is no direct correspon¬ 
dence between the values we employed in our finite- 
element calculations and the parameters used in our LB 
simulations, since we are interested in fitting our simu¬ 
lations to the behavior expected for the Stokes’ flow so¬ 
lution. We used reduced units throughout for our com¬ 
parison between the raspberry particle simulations and 
the COMSOL calculations, therefore any specific choice 
of the parameters is factored out. 


C. Notations Used throughout this Manuscript 

In this section, we summarize the notations used 
throughout this manuscript. This will aid in the un¬ 
derstanding of our results, as many of the notations are 
necessarily similar. 

• H , the plate separation for the confined channel 
system, see Fig. [2ja). 

• Rj/, the effective hydrodynamic radius obtained by 
extrapolating translational mobility measurements, 
see Figs. [2](a,b) , for the limit of plate separation 
H f oo. The subscript h is used to differentiate Rh 
from the raspberry radius R. 

• R r , the effective hydrodynamic radius obtained by 
extrapolating rotational mobility measurements, 
see Fig. [2jc) , for the limit of plate separation 
H f oo. 

• 2 , the vertical position of the sphere with respect 
to the center of the channel (2 = 0), see Fig. H a )- 

• ~max> the maximum distance that a sphere can 
move from the center (z = 0) before making hy¬ 
drodynamic contact with a wall: 2 max = H/2 — R r 
and 2 max = H/2 — R r , respectively, see Fig. H a )- 

• /Iq , the bulk translational mobility. 

• p R , the bulk rotational mobility. 

• Pz ±> the translational mobility for movement along 
an axis perpendicular to the plates as a function of 
2 for fixed plate separation, see Fig. Ha). 

• the translational mobility for movement along 
an axis parallel to the plates as a function of 2 for 
fixed plate separation, see Fig. Hb). 

• p R ± , the rotational mobility for rotation about an 
axis perpendicular to the plates as a function of 2 
for fixed plate separation, see Fig. He) left. 

• p R u, the rotational mobility for rotation about an 
axis parallel to the plates as a function of 2 for fixed 
plate separation, see Fig. H c ) right. 
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• fijj , the translational mobility for movement 
along an axis perpendicular to the plates as a func¬ 
tion of the plate separation H for a sphere in the 
center of the channel z = 0, see Fig. [2ja). 

• fijj ||, the translational mobility for movement along 
an axis parallel to the plates as a function of the 
plate separation H for a sphere in the center of the 
channel z = 0, see Fig. [2^b). 

• , the rotational mobility for rotation about an 
axis perpendicular to the plates as a function of the 
plate separation H for a sphere in the center of the 
channel z = 0, see Fig. [2jc) left. 

• ||, the rotational mobility for rotation about an 
axis parallel to the plates as a function of the plate 
separation H for a sphere in the center of the chan¬ 
nel z = 0, see Fig. |2jc) right. 

• /, the fractional deviation between two results. 

One final point of possible confusion is the use of the 
terms ‘parallel’ and ‘perpendicular’ to the plates. By 
parallel, we mean motion in a direction that is along the 
plates (orthogonal to the normal that defines the plates). 
By perpendicular motion, we mean in the direction par¬ 
allel to this normal, i.e., away from one plate and towards 
the other. For rotation experiments these terms are used 
to describe the orientation of the axis of rotation with re¬ 
spect to the plates. Please refer to Fig. [2] and the above 


list if there is any confusion in reading the text. 

III. RESULTS 

In this section, we discuss the results that we obtained 
by performing the simulations and numerical calculations 
outlined in Section eh 

A. Perpendicular Translational Motion between 
Two Plates 

Let us begin by considering motion perpendicular to 
the plates, see Fig. |2ja). This is the simplest case, since 
there are no translation-rotation cross-coupling terms in 
the HMT for this direction. [35] The results shown in this 
section are for a filled raspberry model with radius R = 
3.0cr and N tot = 925 coupling points. However, we have 
also considered several other radii and the corresponding 
hollow equivalents of these raspberries. We will comment 
on this at the end of the section. 

There are analytic expressions for the positional depen¬ 
dence z of the perpendicular translational mobility 
to which we can compare our simulation results. These 
expressions follow from the result derived by Brenner for 
perpendicular movement near a single wall |9| and were 
extended in Refs. mm to the geometry of two parallel 
plates. N.B. Both Refs. [TD]|TT] contain a typo in formu¬ 
lating the expression in Ref. [9j, the correct expressions 
are given by 


a = cosh 1 (p/R ); 


Ai(p) = -sinh(a)^ 


i(n + 1) 


2 sinh((2n + l)a) + (2 n + 1) sinh(2a) 
^ V(2 n — l)(2n + 3 )) V4sinh((n + l/2)a) 2 — (2 n + l) 2 sinh(a) : 


- 1 


(1) 

( 2 ) 


X 2 (z) = 1 + 5Z 


n =0 L 


Ai ( —(2 n + 1) + z ) — 1 


£ 

n —0 L 


aClO) = A x \p); 

/U,_L = A 2 (~)- 


Ai ( f-(2n-l)-z) -1 


-2^[A 1 (nff)-l]; (3) 


n—0 


(4) 

(5) 


Here, ± {p) is the translational mobility of a sphere 
of radius R moving perpendicular to a single wall. The 
variable p is the minimal distance between the center of 
the sphere and this wall; the parameter a(p) is introduced 
to ease the notation. The perpendicular translational 
mobility between two parallel plates is obtained using 
p-Jj_(p) via the image formalism. N.B. There was a 
missing factor of (4/3) sinh(a) in Eq. (|2j) in our J. 
Chem. Phys. paper, which has been corrected 
here. 

Figure [3](a) shows the result of our quiescent velocity 
and thermalized MSD experiments, Figs. [2|a) and !<*)> 


respectively, compared to the analytic value of p,Tj _, see 
Eq. (§. Here, we used a channel height of H = 14.0(7, in 
which we placed a filled raspberry with radius R = 3.0(7. 
As can be seen in Fig. |3](b) the agreement between our 
quiescent simulations and the corrected series expression 
is quite excellent. There is a fractional deviation / of less 
than 1% throughout the channel, up to the point where 
the raspberry and the channel walls are separated by less 
than 1.0(7. Here it should be noted that the mapping of 
the position to z/z mAx is based on R%. We will come 
back to the effective hydrodynamic radius R^, by which 
we calculated z max , and the bulk mobility by which we 
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FIG. 3. (color online) The translational mobility perpendicu¬ 
lar to the channel walls Zz ± for a sphere of radius R = 3.Oct 
moving through a channel of height H = 14. Oct. The mobil¬ 
ity is normalized by the bulk translational mobility Zo and is 
given as a function of the position in the channel z, which itself 
is normalized by z ma x- The parameter z m ax is the maximum 
distance that the sphere can move from the center (z = 0) 
before making contact with a wall, see the main text for fur¬ 
ther details, (a) The blue circles show the results of quiescent 
LB simulations, the green squares with error bars the results 
of thennalized LB simulations (the error bars indicate the 
standard error), and the dashed red curve the analytic result 
(Eq. © and Refs. mm)- The dashed vertical line (thick 
black) indicates the position for which the sphere and the wall 
are separated by one LB lattice spacing, (b) The fractional 
deviation / of the quiescent LB result from the theoretical 
result. The horizontal gray line shows a fractional deviation 
of 2.5%. 


normalize /jZ ± shortly. 

For small separations there is a systematic deviation 
between the theory and our LB results, see Fig. §a). 
This deviation can be attributed to the break-down of the 
description in terms of an effective hydrodynamic radius 
in proximity of a wall (within one lattice spacing). More¬ 
over, the lubrication forces are not accurately accounted 


for in our combined MD-LB model for close wall-particle 
proximity. This could explain why fJ,Z ^ Z 0 for z = z max . 
Finally, note that the results of our thermalized MSD ex¬ 
periment match the trend of the theoretical curve as well, 
within the error bar. 




FIG. 4. (color online) The translational mobility perpendicu¬ 
lar to the channel walls Zh,± f° r a sphere of radius R = 3.0a 
moving through a channel of height H in the center of the 
channel 2 = 0. The mobility is normalized by the bulk trans¬ 
lational mobility Zo an d is given as a function of the recipro¬ 
cal of the channel height H, multiplied by the hydrodynamic 
diameter (twice the hydrodynamic radius R jC). (a) The blue 
circles show the results of quiescent LB simulations and the 
dashed red curve the semi-analytic result. [91 4111 The dashed 
vertical line (thick black) indicates the position for which the 
sphere (which is in the center of the channel) and the two walls 
are separated by one LB lattice spacing, (b) The fractional 
deviation / of the quiescent LB data from the theoretical re¬ 
sult. The horizontal gray line shows a fractional deviation of 
2.5%. 

In order to determine the value of the effective hy¬ 
drodynamic radius, we varied the distance between the 
plates and measured at z = 0 (the middle of the 
channel), where the mobility attains its highest value. 
We denote this plate-separation dependent mobility by 
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j_. The result of our simulations, compared to the 
theoretical expression in Eq. ([5]) (setting z = 0 and vary¬ 
ing H) is shown in Fig. [4j Agreement between the two 
data sets is quite excellent for H > 2 R% + a, as can 
be appreciated from Fig. |4](b), which shows systematic 
deviation between the theory and numerical results, but 
it is small. The structure in / originates from the fit¬ 
ting procedure used to extract the bulk mobility and its 
subsequent use as a normalization factor. Furthermore, 
the model breaks down for smaller separations, as can be 
clearly seen in Fig. B a )> since sub-lattice changes in the 
plate separation cannot be obtained using the bounce- 
back boundaries. 

Finally, repeating the same experiments with the hol¬ 
low model reveals that the effective radii that we ob¬ 
tained from both models were close to the ones that we 
obtained from our experiments in the simple-cubic ge¬ 
ometry, which we discussed in Part I. (3[ This is shown 
in Table |T] and will be discussed in detail in Section El 
We also performed plate separation simulations for our 
raspberries with R = 2.0ct, 2.5ct, 3.0cr, 4 .Oct, and 5 .Oct. 
For both the hollow and the filled raspberry model and 
all sizes, the dependence of ^ ± and fijj ± on 2 and H , 
respectively, is the same within the error bar (not shown 
here). 


B. Parallel Translational Motion between Two 
Plates 


We continued our investigation into the quality of 
the raspberry model under confinement by examining 
the mobility in the direction parallel to the plates, see 
Fig- lUb). That is, motion perpendicular to the normal 
vector that defines the plate orientation. For such mo¬ 
tion there are translation-rotation cross-coupling terms 
in the HMT. [35] Only when the sphere moves exactly 
in the plane centered between the two plates, there is no 
coupling due to symmetry. 

Figure [5]( a) shows the result of our quiescent velocity 
(Fig. [2](b)) and thermalized MSD experiments (Fig.[2]]d)) 
compared to the /iW obtained using COMSOL. The 
COMSOL results were fitted with a fourth order poly¬ 
nomial. Here, we used a channel height of H = 14.Oct in 
which we placed a filled raspberry with radius R = 3.Oct. 
As can be seen in Fig. [5](c) the agreement between our 
quiescent simulations and the finite-element calculations 
is quite excellent, despite the fact that we ignored cross- 
coupling in the latter. The fractional deviation / is less 
than 1% throughout the channel, even for points where 
the raspberry and the channel walls are separated by less 
than 1.0 ct. The results of our thermalized MSD experi¬ 
ment match the trend of the finite-element curve as well, 
within the error bar, see Fig. [5|a). We should note that 
in our simulations we found no evidence of translation- 
rotation cross-coupling. This is presumably due to the 
small particle-wall separations, for which these effects 
typically express themselves. [35] 
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FIG. 5. (color online) The translational mobility parallel to 
the channel walls ^ y for a sphere (R = 3.Oct) moving through 
a channel of height H = 14.Oct, plotted as a function of the 
position z. The notation is otherwise the same as in Fig. [3] 
(a) Quiescent LB simulations (blue circles), thermalized LB 
simulations (green squares with error bars), and COMSOL 
results (dashed red curve), (b) The fractional deviation / 
between the quiescent LB and COMSOL result (gray line / = 
2.5%). 


For motion in the center plane there is an analytic 
expression for height dependence of the parallel trans¬ 
lational mobility [ijj y derived by Faxen [3 |8j. Faxen’s 
law for center-plane parallel motion between two plates 
is given by 



1.000 - 1.004a: + 0.418a; 3 + 0.210a; 4 - 0.169a; 5 ; 



( 6 ) 

(7) 


It is claimed that Eq. ([6]) is accurate up to fifth order in 
the plate separation. [8] However, the fifth order coeffi¬ 
cient was only provided in a private communication [E] 











and no account of its derivation exist. 


(a) 


I— o 

hi 

zi 


1 

0.8 

0.6 

0.4 

0.2 

0 



- Raspberry o 


Comsol 

Theory 


0 0.2 0.4 0.6 0.8 1 


2R^/H 


(b) 



FIG. 6. (color online) The translational mobility parallel 
to the channel walls in the center of the channel y for 
a sphere of radius R = 3.0<r moving through a channel of 
height H. The notation is otherwise the same as in Fig. [4] 
(a) The blue circles show the data from quiescent LB simula¬ 
tions, the dashed green curve the theoretical approximation 
by Faxen (Eq. (|6|; Refs. mm), and the dashed red curve the 
result of finite-element calculations performed using COM¬ 
SOL (Eq. |8|). (b) The fractional deviation / of the quiescent 
LB result from the theoretical result by Faxen (green squares) 
and the COMSOL result (red circles). The horizontal gray 
line shows a fractional deviation of 2.5%. 


To determine the bulk mobility and the associated ef¬ 
fective hydrodynamic radius, as well as examine the qual¬ 
ity of Faxen’s law, we examined /r J y for parallel mo¬ 
tion in the center plane. The result of our simulations 
is shown in Fig. [6j which provides a comparison between 
our data, the theoretical result by Faxen (Eq. ©), mi 
and the COMSOL finite-element calculations. The re¬ 
sult obtained using the raspberry model agrees perfectly 
with the finite-element data over the entire range of plate 
separations considered here. As expected there is a sig¬ 
nificant deviation between Faxen’s result for small plate 


separations, since the theory ignores higher-order contri¬ 
butions in H to the mobility. We only find acceptable 
agreement between our data and Faxen’s expression for 

H >±R T h - 

From our finite-element calculations we can establish 
an ‘improved’ version of Faxen’s law. The COMSOL data 
points in Fig. |§ a ) were fitted using a fourth-order poly¬ 
nomial 


—= A + Bx + Cx 2 + Dx 3 + Ex a \ (8) 
Mo 

A = +0.997 ± 0.001; (9) 

B = -0.949 ± 0.009; (10) 

C = -0.215 ± 0.034; (11) 

D = +0.895 ±0.051; (12) 

E = —0.542 ± 0.025. (13) 


The numbers with error bars give our numerically ob¬ 
tained coefficients. N.B. The coefficients here are 
correct, for the ones listed in our J. Chem. Phys. 
publication the coefficients D and E are reversed. 

Note that we used a different form of polynomial expan¬ 
sion for our improved Faxen’s law (leaving in the x 2 term) 
to achieve a compacter expression and better \ 2 values, 
since we fit lower-order terms. The first coefficient A can 
be taken to be 1 without significantly changing the result 
of the fitting procedure. This would give the correct limit 
for x 4- 0, which corresponds to infinite plate separation. 

We again found for both the hollow and the filled rasp¬ 
berry, as well as all raspberry radii which we investigated, 
that the dependence of /xj y and /x|] y on z and H , respec¬ 
tively, is the same within the error bar (not shown here). 
The value of the effective radii for the hollow and filled 
models differed, though, see Table [T] This is in agreement 
with our findings for the perpendicular motion between 
plates and in the simple-cubic geometry. [3] We will dis¬ 
cuss this further in Section [TV] 


C. Rotational Motion between Two Plates 

Finally, we considered the accuracy of the raspberry 
model under confinement for the case of rotation about 
an axis in the direction parallel and perpendicular to the 
plates, see Fig. |2](c) . For rotation about an axis paral¬ 
lel to the wall (perpendicular to the normal, see Fig.[2](c) 
(left)) there are rotation-translation cross-coupling terms 
in the HMT. [35] Only when z = 0, there is no such cou¬ 
pling due to symmetry considerations. To the best of our 
knowledge no analytic expressions for rotational mobili¬ 
ties between two plates exist. Therefore, we compare to 
finite-element calculations, the data points of which we 
fitted using a fourth-order polynomial. We only show re¬ 
sults for the filled raspberry model for the position z and 
height H dependence of the mobilities. 

Figure [7](a) shows the result of our quiescent torque 
experiments compared to the data obtained using COM¬ 
SOL for both forms of rotation. Here, we used a channel 
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FIG. 7. (color online) The rotational mobility for rotation 
about an axis perpendicular and parallel to the walls 
and respectively) for a sphere of radius R = 3.Oct mov¬ 
ing through a channel of height H = 14.Oct. The mobilities 
are normalized by the bulk rotational mobility no and plot¬ 
ted as a function of the position 2 from the center (2 = 0), 
expressed in terms of the maximal deviation 2 m ax- The no¬ 
tation is otherwise the same as in Figs. [3] and [5] (a) The 
blue circles show and the blue squares show for 
quiescent LB simulations, the dashed red and green curves 
show the corresponding results obtained using COMSOL. (b) 
The fractional deviation / between the quiescent LB and the 
COMSOL result, rotation about an axis perpendicular to the 
walls is indicated using red circles and rotation about an axis 
parallel to the wall using green squares (gray line / = 2.5%). 


height of H = 14.0er in which we placed a filled rasp¬ 
berry with radius R = 3.0er. As can be seen in Fig. m 
the agreement between our quiescent simulations and the 
finite-element calculations is quite excellent. The frac¬ 
tional deviation / is less than 1% throughout the chan¬ 
nel for rotation about an axis parallel to the wall and 
even better for rotation about an axis perpendicular to 
the wall. We can see significant deviation for the former 
when the raspberry and the channel walls are separated 


by less than l.Oer, but this can be explained by the small 
particle-wall separation compared to the lattice spacing. 
Again, we did not observe translation-rotation cross cou¬ 
pling in our data. 

Note that there is some structure to the deviation be¬ 
tween the parallel-rotation raspberry data and the finite- 
element calculations. This is related to the fitting pro¬ 
cedure used to extract the bulk rotational mobility q 
by varying H. A small deviation in the effective hydro- 
dynamic radius Rff can have significant impact on the 
shape of the curve, as we will explain shortly. For the 
rotation about an axis perpendicular to the wall there is 
a similar increase when the wall-raspberry separation is 
less than 1.0<r, but the agreement is still quite reasonable. 

To determine the bulk mobility and the associated ef¬ 
fective hydrodynamic radius, we examined /zjj j_ and n 
(rotation in the mid-plane with varying H), see Fig [8] 
Again we find quite acceptable agreement between the 
data obtained using the raspberry model and the finite- 
element calculations. However, we can clearly see that 
there is structure to the fractional deviation / between 
the two results in Fig. |s](c) . This can be explained as 
follows. 

In determining the effective hydrodynamic radius, we 
used a quartic polynomial to fit our data, which was given 
as a function of 2R^/H and normalized by the value of 
Hq corresponding to that R That is, the data to which 
we fitted is a function of R and we applied a root-finding 
algorithm to ensure that the fitted value for H f oo of 
Vh,.l/Ro an d Rh \\/Ro equals 1. We thus solved for the 


value Rfr- 

The above procedure, however, strongly favors the cri¬ 
terion imposed on the extrapolated data point, while it 
sacrifices overall matching to the finite-element calcula¬ 
tions - which we assume to be more accurate. Hence, 
there can be a systematic deviation between the two data 
sets, as is indeed found in Fig. |8jc). However, our fitting 
procedure is reasonable, since it does not implicitly as¬ 
sumes the finite-element data to be superior to that ob¬ 
tained using the raspberry model and it only uses internal 
fitting parameters. The agreement between the two re¬ 
sults, despite the systematic deviations, can therefore be 
considered excellent. 

From our numerical results we can formulate modified 
empirical forms for the rotational mobility’s dependence 
on the channel height H , for mid-channel rotation about 
an axis parallel and perpendicular to the walls. These 
expressions are as follows 


Rha 

Ro 


= 1 + AiX 2 + BiX 3 ± CiX 4 


(14) 


A± 
B± 
C _L 


2 Rl. 

H ’ 

- 0.10 ± 0 . 01 ; 
±0.09 ±0.03; 
-0.26 ±0.02; 


Ay = -0.09 
H|| = -0.17 
Cy = -0.38 



(15) 

±0.01; 

(16) 

± 0.04; 

(17) 

±0.03, 

(18) 


where i is a subscript that can be either ± for rotation 
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FIG. 8. (color online) The perpendicular (a) and parallel (b) 
rotational mobility (/x:? x and /x?||i respectively) for a sphere 
of radius R = 3.Oct moving through a channel of height H. 
The mobilities are normalized by the bulk coefficient /xj and 
plotted as a function of 1 /H, expressed in terms of twice the 
hydrodynamic radius R ^ . The notation is otherwise the same 
as in Figs. [4] and [6] Quiescent LB simulations are given by 
blue circles and COMSOL results using dashed red curves, (c) 
The fractional deviation / of the quiescent LB results from 
the COMSOL data (gray line / = 2.5%). 


about an axis perpendicular to the walls and || for rota¬ 
tion about an axis parallel to the walls. We divided out 


the linear coefficient /x^ to obtain the functional form of 
our result. 

Again, we found that for both the hollow and the filled 
raspberry model, as well as all three raspberry radii which 
we studied, the dependence of /x^m, /x^| j_, and /x^ n 
on z and H, respectively, is the same within the error 
bar (not shown here). The value of the effective radii 
for the hollow and filled models differed in each case, see 
Table[IJ in agreement with our findings for the perpendic¬ 
ular and parallel motion between plates, as well as in the 
simple-cubic geometry. )3j We will discuss this further in 
Section llVl 


IV. DISCUSSION 

In Section m we have shown that there is excellent 
agreement between established analytical and numerical 
results for the Stokesian hydrodynamic behavior of spher¬ 
ical particles compared to their raspberry equivalent in 
the confining geometry of two parallel plates. This ex¬ 
tends our findings of Part I [5J to geometries with walls. 
The discrepancy between the translational and rotational 
hydrodynamic radius of the original hollow raspberry - 
first observed by Ollila et al. [5] - is present for the par¬ 
allel plate system as well. In this section we discuss this 
discrepancy in more detail. 

In Table |T] we summarize the results for the hy¬ 
drodynamic radii obtained in our experiments in this 
manuscript as well as in Part I. [3] Note that within 
the error bar the data for the filled raspberries are fully 
consistent, whereas the data from the hollow raspber¬ 
ries show a clear discrepancy between R ^ and R The 
Rh obtained from the experiment (R = 3.Oct) with rota¬ 
tion about an axis perpendicular to the wall (Fig. [2jc) , 
right) seems to be an outlier for both the filled and hol¬ 
low raspberry, but it is still within two standard devia¬ 
tions of the mean for the filled raspberry. As explained 
in Section |III C[ our fitting procedure ensured that for 
the box length it 00 ' 11 the simple-cubic geometry 0 
and for the plate separation H f oo in the confined ge¬ 
ometry, the fitted result converges to the expected bulk 
mobility for a given hydrodynamic radius, for which we 
subsequently solved. Therefore, the rotational data has 
a higher standard deviation in general, as taking the cu¬ 
bic root to establish the effective radius leads to stronger 
error propagation. 

From our results, we conclude that the level of dis¬ 
agreement between the effective radii of the hollow rasp¬ 
berry model seems to be unaffected by the confinement. 
This can be seen from the excellent agreement (in an ab¬ 
solute sense) between the hydrodynamic radii in bulk and 
good relative agreement between our results and litera¬ 
ture references in the confining geometry. We have also 
shown that our ‘filling + fitting’ procedure for the rasp¬ 
berry model (see Ref. 0) again significantly improves 
the agreement between the effective hydrodynamic radius 
obtained by translational and rotational experiments. 
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filled 

hollow 

R = 2.0cr 

Rl 

Rl 

Rl 

pi? 

simple-cubic 

2.59cr 

2.57 a 

2.50 a 

2.42a 

walls ± 

2.57a 

2.56a 

2.48a 

2.43a 

walls | 

2.57a 

2.56a 

2.48a 

2.43a 

R = 2.5cr 

Rl 

p i? 

Rl 

pi? 

simple-cubic 

3.01cr 

3.02 a 

2.98 a 

2.90a 

walls ± 

3.00cr 

3.01a 

2.97a 

2.90a 

walls || 

3.00cr 

3.01a 

2.97 a 

2.90a 

R = 3.0(7 

Rl 

pi? 

Mil 

Rl 

pi? 

simple-cubic 

3.53cr 

3.54a 

3.47a 

3.38a 

walls ± 

3.53cr 

3.50a 

3.48a 

3.44a 

walls || 

3.53a 

3.51a 

3.47a 

3.37a 

R = 4.0cr 

Rl 

pi? 

Rl 

pi? 

simple-cubic 

4.57a 

4.56a 

4.46a 

4.37a 

walls ± 

4.55a 

4.54a 

4.45a 

4.36a 

walls || 

4.54a 

4.54a 

4.44a 

4.36a 

R = 5.0cr 

Rl 

Rl 

Rl 

pi? 

simple-cubic 

5.55 a 

5.53 a 

5.44a 

5.35 a 

walls ± 

5.55a 

5.52a 

5.44a 

5.35a 

walls || 

5.54a 

5.52a 

5.43a 

5.35a 


TABLE I. Summary of the various effective radii determined 
in the experiments described in Section III and Ref. [JJ. There 
is a ±0.01 ct error bar for the translational experiments and 
±0.02cr for the rotational experiments. The bare friction co¬ 
efficient used to generate this data is £o = 25moT _1 . 


There is correspondence between the effective radii that 
are found for the ‘filling ± fitting’ raspberry using the 
various experiments, within the numerical error. A simi¬ 
lar conclusion cannot be reached for the hollow raspberry, 
indicating that this model may be more sensitive to the 
details of its surroundings. 

Finally, let us focus on the raspberry constructed us¬ 
ing the Ahlrichs and Diinweg viscous coupling (VC) 
scheme [5] in the context of the other methods avail¬ 
able to model the hydrodynamic interactions of col¬ 
loids with each other and boundaries. There are sev¬ 
eral well-known alternatives to the method outlined in 
Part I [3] and this manuscript, by which MD objects can 
be coupled to LB, namely: Ladd bounce-back bound¬ 
aries (Ladd BB), [M1EZ] the immersed boundary method 
(IBM), [55] and the external boundary force (EBF) 
method. J5S] There are other hydrodynamics solvers that 
can achieve similar colloid-in-fluid descriptions. Most 
commonly used in physics are dissipative particle dy¬ 
namics (DPD), [10] jH] multi-particle collision dynamics 
(MPCD) or stochastic rotation dynamics (SRD), [421143] 
and Stokesian Dynamics (SD). [44] In the following we 
briefly discuss our method in the context of these alter¬ 
natives. 

It was recently shown by Schiller [35] that VC is a 
general description for fluid-particle coupling, which en¬ 


compasses both IBM and EBF. IBM and EBF corre¬ 
spond to special choices of the friction and mass ratio 
in VC. As such, for the LB-based descriptions it suffices 
to consider Ladd BB. One point of interest is that the 
EBF method can be modified to lead to an instantaneous 
no-slip boundary condition, but only at the expense of 
loosing freedom in the choice of the bare friction coeffi¬ 
cient. [35] 

Ladd BB employs the grid on which the LB fluid is 
simulated to describe the particles. [Ml EH This gives 
Ladd BB an advantage over the VC raspberry model, 
as it does not necessitate an interpolation step or sub¬ 
lattice refinement to achieve the fluid-particle coupling. 
Similar to the VC raspberry particle, Ladd BB also gives 
rise to an effective hydrodynamic radius, for which can 
be fitted. However, Ladd BB achieves a no-slip boundary 
condition, without necessitating a particular choice of the 
bare friction coefficient. This is favorable, as there is no 
question of porosity. [5H5] It is therefore likely that Ladd 
BB performs better in force and velocity experiments of 
Part I, [3] as no counter-force can be applied to the nodes 
which are encompassed by the Ladd BB particle, leading 
to an effective (reduced) force being applied to the par¬ 
ticle. Finally, there are lubrication corrections available 
for the Ladd BB method. [Ml 07] These should improve 
the near-field interactions between particles and particles 
and walls, compared to the VC raspberry method. 

Considering these advantages Ladd BB would seem the 
method of choice. However, the grid-based description of 
the particle in Ladd BB is known to suffer when there 
are too few grid points inside the particle. For exam¬ 
ple, to obtain a coupling shape that is sufficiently spher¬ 
ical, a radius of R > 2.5(7 is required. We have shown 
that the VC raspberry can be used for smaller R and 
the VC has been shown to work well even for single cou¬ 
pling points. [2j This is useful in systems where one is 
restricted in the number of LB nodes that one can use 
- in particular graphics processor unit (GPU)-based LB 
codes are limited by the amount of RAM available on the 
GPU. VC allowes one to study large systems with many 
thousands of particles, without the need to fully resolve 
each on the lattice, see, e.g ., Ref. [351 . 

The non-LB methods to achieve fluid-colloid cou¬ 
pling vary considerably in the way that this coupling is 
achieved. DPD, MPCD, and SRD (the latter two essen¬ 
tially being the same), all achieve coupling to the ‘fluid’ 
by interacting with the particles that comprise it. For 
MPCD and SRD, the smallest length scale on which the 
larger raspberry object can be resolved is governed by 
the size of the fluid particles, rather than the boxes that 
are used in the collision step. This could well be advanta¬ 
geous in modeling small separations of walls and raspber¬ 
ries. However, it should be noted that a sufficient density 
of MPCD particles must be present to properly model 
hydrodynamic interactions, which would be problematic 
for situations in which the raspberries are in proximity 
to each other and obstacles. In DPD, however, the finest 
features of the raspberry that can be resolved have the 
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typical length scale of the solvent particle interactions. 
As such, DPD is typically a coarser method than both 
the VC in LB and MPCD and is more likely to encounter 
problems for small raspberry-raspberry and raspberry- 
wall separations. 

Finally, there is the method of SD, [H] in which the hy¬ 
drodynamic interactions between particles are described 
using the Rotne-Prager tensor. This method a suited for 
bulk systems, but it can be modified to periodic bound¬ 
ary conditions. Its main attraction is the fact that it 
solves Stokes’ equations exactly for dilute systems and 
that it can be modified to include a lubrication correc¬ 
tion. This together with the algorithm’s favorable scal¬ 
ing [IS] give it a slight edge over the VC method for 
bulk systems. However, simple (closed) analytic forms 
for the Rotne-Prager tensor are not available for more 
complicated geometries - not even for two plates. Only 
for single plates can one use the Blake-tensor formalism. 
Therefore, any advantage of the SD algorithm in bulk 
(with or without lubrication corrections) is lost in more 
complicated geometries, since the quality of the solution 
sensitively depends on the approximations made for the 
hydrodynamic interaction tensors. 

A full comparison of the various methods goes beyond 
the scope of the current research, but is certainly worth¬ 
while for future study. 


V. CONCLUSION AND OUTLOOK 

Summarizing, we have examined the properties of the 
raspberry model using a classic fluid dynamics experi¬ 
ments, wherein a sphere is subjected to the confining 
geometry of parallel plates. This so-called ‘raspberry’ 
model refers to a hybrid lattice-Boltzmann (LB) and 
Langevin molecular dynamics (MD) scheme for simulat¬ 
ing the dynamics of suspensions of colloids originally de¬ 
veloped by Lobaskin and Dunweg. [T] The particle is rep¬ 
resented by a set of points on its surface that couple to 
the fluid through a frictional force acting both on the 
solvent and on the solute, which depends on the relative 
velocity. 

We began our investigation by analyzing the position 
dependent translational mobility of a sphere moving per¬ 
pendicular to the two plates. We found excellent agree¬ 
ment with the predictions of Refs. mm. This was ex¬ 
tended upon by considering the parallel translational mo¬ 
bility of the sphere. Here, we found good agreement with 
the result by Faxen [313 in the limit of large plate sep¬ 
aration. We improved upon Faxen’s result by propos¬ 
ing a modified, empirical expression that works well over 
the entire range of plate separations. Finally, we con¬ 
sidered the rotational mobility for rotation about axes 
parallel and perpendicular to the plates. For these sys¬ 
tems no analytic expressions were available to which we 
could compare our numerical results. However, on the 
basis of our numerical results, we were able to formu¬ 
late expressions for these mobility coefficients for spheres 


fixed in the center of the plates as a function of the plate 
separation. 

Our results show that Stokesian hydrodynamic behav¬ 
ior is reproduced by our raspberry model with the ‘filling 
+ fitting’ formalism to a surprising degree of accuracy 
over a wide range of length scales. From our combined 
data we can draw the following additional conclusions 
concerning our ‘filling + fitting’ procedure and the rasp¬ 
berry point-coupling method in general: 

• To determine the bulk translational mobility of a 
particle, a force or velocity experiment performed 
in the center of two plates is more suitable than 
in a simple-cubic geometry. This type of exper¬ 
iment does not have the disadvantage of a back 
force/velocity density, [3 which could interfere with 
the fitting procedure required to extrapolate results 
to bulk (infinite plate separation). 

• The raspberry point-coupling can be combined 
with confining geometries consisting of bounce- 
back boundary conditions. The combination yields 
accurate results over a large range of particle¬ 
boundary separations, up to roughly one LB lattice 
spacing; as verified for the specific geometry of two 
parallel plates. 

• We estimate the raspberry to accurately reproduce 
hydrodynamics interactions for particle-particle 
and particle-surface separations greater than one 
LB lattice spacing. The separation is determined 
by the effective radii. Here, we consider the results 
of our rotational mobility experiments and transla¬ 
tional mobility experiments in confinement signif¬ 
icant, as there are no back-torques or back-forces 
applied to the fluid. [5] 

The latter estimate is not unreasonable, since the near¬ 
field effects of discretizing the surface into a finite number 
of MD beads become noticeable at comparable separa¬ 
tions. There are, however, two caveats to the last result, 
(i) We have not verified that we obtain the correct be¬ 
havior in a system with many particles that can undergo 
collision events. Here, we only extrapolate our particle- 
image results 0] to particle-particle interactions. This 
point is left for future study, (ii) In the case of a rotat¬ 
ing raspberry, for which we observed that the particle- 
image interactions are correct up to one lattice spacing, 
the motion and structure of the raspberry are such that 
the fluid is constantly exposed to different fluid-coupling 
points (the MD beads). This may remove artifacts caused 
by the discretization of the surface (or volume) that are 
more noticeable in the translational experiment, where 
the orientation is fixed. 

From the above and our examination in Part I [5], it 
becomes clear that our ‘filling + fitting’ procedure for 
the raspberry model is an excellent way to approximate 
fluid-particle coupling in a Stokes’ liquid, using an LB 
algorithm. 
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